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ABSTRACT 



We present the results of near- to mid- infrared slit spectroscopic observa- 
tions (2.55-13.4/im) of the diffuse emission toward nine positions in the Large 
Magellanic Cloud with the Infrared Camera (IRC) on board AKARI. The tar- 
get positions are selected to cover a wide range of the intensity of the incident 
radiation field. The unidentified infrared bands at 3.3, 6.2, 7.7, 8.6 and 11.3/im 
are detected toward all the targets, and ionized gas signatures: hydrogen re- 
combination lines and ionic forbidden lines toward three of them. We classify 
the targets into two groups: those without the ionized gas signatures (Group A) 
and those with the ionized signatures (Group B). Group A includes molecular 
clouds and photo-dissociation regions, whereas Group B consists of H II regions. 
In Group A, the band ratios of /3.3/i,m//ii.3/.m, h.2i,m/hi.3f,m, h.r i^m/ hi.si^m and 
-^8.6 /im /-^11.3/im show positivc Correlation with the IRAS and ^i^^Ai?/ colors, but 
those of Group B do not follow the correlation. We discuss the results in terms 
of the polycyclic aromatic hydrocarbon (PAH) model and attribute the differ- 
ence to the destruction of small PAHs and an increase in the recombination due 
to the high electron density in Group B. In the present study, the 3.3 fim band 
provides crucial information on the size distribution and/or the excitation con- 
ditions of PAHs and plays a key role in the distinction of Group A from B. The 
results suggest the possibility of the diagram of h.Sf^m/ hi.3f,m v.s. I7.7 f,m/ lu.Sf^m 
as an efficient diagnostic tool to infer the physical conditions of the interstellar 
medium. 

Subject headings: dust, extinction — Magellanic Clouds — galaxies: ISM — infrared: 
galaxies — infrared: ISM 
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Introduction 



Since the discovery of tlie 11.3yum band in planetary nebulae in 1973 flGillett et al. 



19731 ). space infrared missions, and air-borne and ground-based infrared observations have 
shown that the major unidentified infrared (UIR) bands appear at 3.3, 6.2, 7.7, 8.6, 11.3, 
12.6 and 16.4 fim together with some faint features. The UIR bands have been observed in 
various astrophysical environment s, including photo-d issociation regions (PDRs), reflection 



nebulae, pla netary nebulae (e. g.. 



(ISM) 



(e-g-. 



Onaka et al. 



Dale et al 



Lutz et al 



2005 



2006 



1996 



Peeters et a. 



Mattila et al. 



Kaneda et al. 



Sajina et al. 



20081 



20021) . the diffuse interstellar medium 



19961). nearby galaxies of various types 



Smith et al. 



20071 ). and distant galaxies (e.g. 



20071 ). The carriers of the UIR bands are generally thought to 



be po^ 



ycyclic arornatic hydrocarbons (PAHs) or PAH-containing car 



Sakata et al. 



198 



Puget fc Leger 



1989 



Papoular et al. 



1989 



jonaceous compoun ds 



Allamandola et al. 



19891 ). PAHs are excited by absorbing a single UV photon and emit a number of IR 
photons corresponding to vibration modes of C-C and C-H bonds. The 3.3 /im band is 
assigned to C-H stretching modes, the 6.2 /xm band to C-C stretching modes, the 7.7 /xm 
band to blending of several C-C stretching modes and C-H in-plane bending modes, the 
8.6 /im band to C-H in-plane bending modes, the 11.3/im band to solo C-H out-of-plane 



bending modes, anc 



( lAUamandola et al. 



the 12.6 /im band to trio C-H out-of-plane bending modes, respectively 



19891 ). 



Recent laboratory experiments and quantum chemical calculations suggest that the 
properties of the UIR bands (e.g., shapes, center wavelengths, interband ratios, etc.) reflect 
the chemical and physical properties of PAHs (e.g., molecular structure, size distribution, 
ionization stat e, temperatur e, etc.), which may be altered in interstellar and circumstellar 



environments (Tielens 



20081 ). Therefore the UIR bands have a great potential to be used as 



efficient diagnostic tools to infer the physical condition of the ISM even in remote galaxies. 
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Observations of the diffuse Galactic radiation and normal galaxies have shown very 
l ittle variation in the mid-infrared (MIR) UIR ban d spectra (6-12 fim) until recently 



(jChan et al. 



2001 



Lu et al 



2003 



Sakon et al. 



20041 ). whereas small variations in the MIR 



UIR bands have be en reported between the disk and halo reg ions or the arm and interarm 



regions of galaxies (jlrwin fc Madden 



2006 



Sakon et al 



20071 ). Latest Spitzer and AKARI 



observations c 


early 


(Kaneda et al. 


2007 



early show distinct variations in the M 



Smith et al 



2007 



Galliano et al 



■j UIR spectra in particular galaxies 



20081 ) for the first time. However, 



the diagnostic for the physical conditions of the ISM as well as the chemical and physical 
evolution of PAHs in galaxies by means of the UIR bands is not yet fully explored. 

In this paper, we present the results of near-infrared (NIR) to MIR spectroscopic 
observations of the ISM with different radiation conditions in t he Large Magellanic 



Cloud (LMC) with the infrared camera (IRC) onboard AKARI (IMurakami et al. 



2007 



Onaka et al. 



2007bl ). The L MC is a near 



50 kpc from the Milky Way (IFeast 



1999 



jy irregular galaxy an d is located at the distance of 



the almost face-on orientation (i ^35*^ 



Keller fc Wood 



20061). In addition to its proximit y. 



van der Marel fc Cioni 



2001 



Nikolaev et al. 



Olsen fc Salyk 



2002 



20041 ) provides us with a unique opportunity to investigate regions with 
different physical conditions without confusion because of the spatial resolution of ~ 5" in 
the MIR of AKARI/mC The IRC spectroscopy has a unique characteri stic that it cari 



obtain a spectrum from 2.5 to 13 fim simultaneously with the same slit (lOhyama et al. 
20071 ). This has an advantage for the spectroscopy of extended objects over other space 
instruments since the Short Wavelength Spectrometer (SWS) onbo ard the Infrared Space 



Observatory (ISO) had different diaphragms from the NIR to MIR fide Graauw et al 



and t h e Infrared Spe c trogra ph (IRS) on Spitzer lacks a channel in the NIR (IHouck et al. 



20041 ). 



Vermeij et al. 



19961) 



(120021 ) report the MIR UIR band ratios of H II regions of the LMC 
based on observations with the ISOPHOT/PHT-S instrument on board ISO, but do not 
include the 3.3 /im band in their analysis because of the low signal-to-noise ratio (S/N) in 
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the short w avelength channeL The 3.3 /im band, in fact, is most sensitive to the smallest 



PAHs (e.g., 



Schutte et al 



1993[ ) and its relative intensity to the MIR UIR bands provides us 
with significant information on the average temperature of PAHs, which depends on the size 
distribution and the excitation condition of PAHs. In this paper, we investigate variations 
in the relative intensity of the UIR bands in the NIR to MIR of the diffuse radiation from 
regions with different radiation field conditions in the LMC and discuss them in relation to 
the physical properties of PAHs. 

In §2, the observation and the data reduction are described together with the selection 
of the target positions. The obtained spectra are presented in §3. In §4, the observed 
variations in the UIR band ratios in different radiation field conditions are investigated in 
terms of the PAH model. Diagnostic of the physical conditions of the observed regions is 
also discussed based on the UIR band ratios. A summary and conclusions are given in §5. 



2. Observations and Data Reduction 
2.1. Observations 

The present study employs the datasets of eight pointed observations (observation IDs: 
1400330, 1400346, 1400318, 1402426, 1400324, 1402422, 1400334 and 1400320) collected as 
part of the AKARI mission program "ISM in our Galaxy and Nearby Galaxies" (ISMGN; 



Kaneda et al 



20091 ). All of the observations were performed with th e slit spectroscopic 



Ohyama et al. 



20071). 



mode with the choice of the grism for the NIR disperser (AOT04 b:Ns 
The NIR spectrum was taken with the grism, NG (2.5-5.0 /im, A/AA ~ 100), in the NIR 
channel and the MIR spectra were taken with two grisms, SGI (4.6-9.2 /im, A/AA ~ 50) 
and SG2 (7.2-13.4 /im, A/AA ~ 50), in the MIR-S channel. The NIR and MIR spectra 
were obtained simultaneously by means of the beam splitter, which enabled us to get a 
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continuous NIR to MIR spectr um from 2.5 to 13.4 fim of the same slit area of 1' length by 



5" width (jOhyama et al 



20071 ). The observation parameters are summarized in Tabled! 
The accurate slit position is determined from the 3.2 fim image (N3), which is taken during 
each pointed observation (see §2.3.ip . by referring to the positions of point sources in the 



2MASS catalog. 



2.2. Target Selection 



The targets are selected by taking account of the CO mapping data (IMizuno et al. 



200ll ) and the IRAS colors of /?5,,m//i^ 



formation activities ( iBoulanger et al 



1988 



and J, 



fin ji.m 



11 



Onaka et al 



^n /"Ti 



, w hich indicate local star 



2007ah 



Sakon et al 



mm have 



shown that the extremely large IRAS colors of h^fim/Iufim and Im^lm/Iloo^lm (~3 and 
~0.6, respectively) in the diffuse emission in the LMC can be accounted for by a large 
contribution from nearby young (< 30Myr) clusters to the incident interstellar radiation 



field and that the small IRAS colors of h^nm/ Ii2nm and /eo^tm/-^: 



100 /iin 



A and ~0.3, 



respectively) are those expected from the heating by the inci dent radiation field inside 



quiescent molecular clouds (e.g. 



Miville-Deschenes et al. 



20021). According to these results, 



we select several infrared bright positions with different IRAS colors as the targets o . 



present study, wh ere molecular clouds are recognized on the CO maps flMizuno et al 



Fukui et al. 



2008 



the 



2001 



). Since the beam size of the IRAS data is larger than the size of the slit of 



Ita et al. 



the A KARI I lYlC, we also derive the AKARI color of /l24/-?^sii from the dataset of 
( I2OO8I ). where /l24 is the flux density at the L24 (24 /im) band and Jgn is that at the Sll 
band (11 /xm) of the AKARI/IRC, to obtain the local radiation field conditions. The flux 
density is measured over an aperture of 5" in diameter around the central position of the 
slit. The IRAS and AKARI colors of the present targets are summarized in Table |2j The 



trend of the AKARI color of /l24/-^sii is very similar to that of the IRAS /25/im/A2/im color. 
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whereby we confirm that the selection based on the IRAS colors is in fact relevant to the 
purpose of the present study. Hence, the target positions cover a wide range of incident 
radiation field conditions, including molecular clouds, PDRs and H II regions. 



In Figure [T], the slit positions of those datasets are shown o ver the fa 



se co lor image 



20081). The Sll 



of the LMC obtained by the AKARI IRC LMC survey program (llta et al.l 
band images of a 10' x 10' area including each slit position are also shown in Figures [2^-h. 
The slit is positioned at regions without apparent point sources in all the positions except 
for Position 8. At Position 8, a bright point-like source is recognized in part of the slit. 
We split the spectrum of Position 8 into 2 parts: the one including the point-like source 
(Position 8-1) and the other without the point-like source (Position 8-2). For the other 
positions, the spectrum is extracted over the almost entire slit length of 40-50". 

The .4ii:^i?/ color of h2i/hii and the IRAS colors of I25f,m/h2f,m and /eo/xm/AooMm at 
Positions 1, 2, 3, 4, 5 and 6 exhibit only a limited range of 0.5-1.8, 1.0-1.9, and 0.3-0.5, 
y . Th ese targets are not associated with SWB or I type star clusters in 



Bica et al. 


(] 


.996 


). 


(Sakon et al. 


2006) 



19961). indicating that they are surrounded by relatively quiescent environments 



values. These targets are located in regions associated with OB star clusters as well as 
H erbig Ae/Be star clust ers (N158-01 and N158-Y1 at Position 7 and N159-Y4 at Position 



Nakajima et al. 



20051 ). confirming that the infrared colors manifest recent star-formation 



activities in the ISM. 



2.3. Data Reduction 



The present data reduction basically follows the standard toolkit for the IRC 
spectroscopy. However, most of the present targets are faint and require careful data 
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processing. Thus, some part of the process is carried out separately from the toolkit with 
special care. Details of the data reduction process are described in the following. 



2.3.1. Slit spectroscopy with AKARI IRC/NIR 

During a single pointed observation with AOT04 grism mode, eight to nine exposure 
frames of NIR grism spectroscopic (NG ) data and one expo sure frame of 3.2 fim imaging 



(N3) data are taken with the IRC/NIR (jOhyama et al. 



20071 ) . The dark current is measured 



in one frame each at the first and last parts of the pointed observation. A single exposure 
frame consists of one short-exposure and one long-exposure images. In the present study, 
only the long exposure data are used. The dark image for each NIR observation is obtained 
by averaging three long-exposure images of the dark current, which are collected from the 
adjacent pointed observations including itself to correct for any high-energy ionizing particle 
(hereafter cosmic-ray) effects by a l.S-a-clipping method. In the present analysis, only the 
dark current data measured in the first part of each pointed observation are used to avoid 
the latent image effects and subtracted from the observation images. 

We recognize small shifts in position due to the pointing instability by at most ~ 5" 
during each pointed observation except for Position 5, where an extraordinary large shift 
of ~ 15" is found. The shift in the direction parallel to the slit is corrected so that the 
spectra of the same area of the sky are extracted. Because the shift in the orthogonal 
direction is uncorrectable, the exposure frames shifted in the direction perpendicular to the 
slit by more than a pixel (~ 1.5") as well as those affected by severe artifacts are discarded, 
except for Position 5 (see §2.3.3). The remaining images are averaged taking account of the 
shift in positions in the direction parallel to the slit by a 1.5-(7-clipping method to remove 
cosmic-ray events and the artifacts. 
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The most critical part in the data-reduction of NG spectra is the removal of artificial 
patterns as well as the foreground components originating from the zodiacal light and 
diffuse Galactic emission. In NIR observations, pixels saturated by cosmic-ray hits or 
extremely luminous objects often produce arti ficial line-like patte rns, sometimes termed as 



"column pull-down" and "multiplexer bleed" (jPipher et al 



20041 ) ■ Particularly an artificial 



line-like pattern in the direction parallel to the slit mimics emission or absorption features in 
the slit spectrum. The position of the artificial line-like structure sometimes differs among 
different exposures even in the same pointed observation. When an artificial line structure 
emerges in a certain exposure image, it is removed by replacing the data of the affected 
pixels with those of the unaffected pixels of other exposure images at the same position. 
To estimate the foreground components from the zodiacal and diffuse Galactic emission, 
we employ the datasets obtained at positions off the LMC with AOT04 a;Ns (observation 
IDs: 1500719 and 1500720). In these observations the NIR spectra are taken with the prism 
mode (NP) instead of the grism mode (NG). Because the emission at the off-LMC positions 
are too faint for the observations with the NG mode, we use the data with the NP mode to 
obtain a reliable foreground spectrum. The slit positions of both observations are centered 
at («2ooo! ^2ooo) = (OG'^OO'^OO.'^O, — 66°36'30"), which is almost at the same ecliptic latitude 

{(3 - -90°) as that of the LMC ( 85°), but is at - 9.7° away from the center of the 

LMC. The observation log of the off-LMC position is also given in Table [H 



2.3.2. Slit spectroscopy with AKARI IRC/MIR- S 

The data reduction procedures including the dark-current subtraction and the 
cosmic-ray correction for MIR-S spectroscopic observations are basically the same as those 
for NIR spectroscopic observations. During a single pointed observation with AOT04, four 
exposure frames of SGI data, four to five exposure frames of SG2 data, and one exposure 
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frame of 9.0 /im imaging (S9W) data are taken with the MIR-S channel. The dark current 
is measured in one frame each in the first and the last parts of the pointed observation. A 
single exposure frame consists of one short-exposure image and three long-exposure images. 
The dark image for each MIR-S observation is obtained by averaging three long-exposure 
images of the dark current by a 1.5-cr-clipping method to correct for the cosmic-ray effects. 
In this process, only the dark current data measured in the first part of each pointed 
observation sequence are used to avoid the latent image effects. The same shifts in position 
as those recognized in the NIR data are expected in the MIR-S observations because the 
same field-of-view is shared by the NIR and MIR-S channels and they are corrected in the 
same manner as in the NIR data. The shift in position among three long-exposure images 
taken in a single frame is negligible in most cases and, therefore, they are averaged by a 
1.5-(T-clipping method to remove the cosmic-ray effects. 

The subtraction of the foreground components (zodiacal light and diffuse Galactic 
emission) is a more serious problem in the data reduction of MIR spectroscopy than in NIR. 
In addition, the MIR de tector suffers scat tered light originating in the scattering within 



the detector array (e.g.. 



SakonetaL 



20071 ). To estimate the foreground emission and the 
scattered light component, the SGI and SG2 spectra collected at the position off the LMC 
(Observation IDs: 1500719 and 1500720) are used. The spectrum at off-position is obtained 
by averaging the two observations and is subtracted from the spectra of the target positions. 
The MIR-S spectra are basically dominated by the zodiacal light and thus the scattered 
light component of the off-position spectrum is almost the same as that in the spectra of 
the target. Therefore, the subtraction of the off-position spectrum works effectively not 
just to remove the foreground emission but also to correct for possible artifacts and greatly 
improves the resultant spectra. 
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2. 3. 3. Continuous spectra from NIR to MIR 

Using the spectral response function of each module, we obtain NG, SGI and SG2 
segmental spectra at the same area of the sky except for positions 5 and 8-1 (see below). 
Each segmental spectrum is truncated at the wavelengths where the S/N becomes low: 
NG is truncated at 2.55 /im and 4.9 /xm, SGI at 5.5 /im and 7.9 /xm and SG2 at 7.9 /xm 



and 13.4 yum. Then, the correction for the slit efficiency for extended sources (jSakon et al. 



20081 ) is applied and continuous spectra from 2.55 to 13.4 /xm are obtained with a small 
gap between 4.9 and 5.5 fim. Because of the severe artifacts (column-pulldown) the NG 
spectrum is truncated at 3.8 fim and 4.5 fim for Positions 2 and 5, respectively. Note that 
although there is a small gap between the NG and SGI segments, all the NG, SGI, and 
SG2 segmental spectra are smoothly connected to each other without scaling, suggesting 
that the subtraction procedure of the foreground emission and scattered light works well 
and rehable spectra are obtained. 

For Positions 5 and 8-1, we cannot obtain segmental spectra precisely at the same 
region of the sky between SGI and SG2 because the SGI and SG2 observations are not 
carried out simultaneously during a pointed observation. A relatively large positional shift 
(~ 15") recognized during the pointed observation for Position 5 prevents us from obtaining 
the SGI and SG2 spectra from the same region of the sky. For Position 8-1, the positional 
stability during the pointed observation is almost the same as the other observations. 
However, a small shift in the position of the point-like source in the slit of 5" width between 
the SGI and SG2 observations changes the source flux to some extent and makes a small 
difference in the flux level between the SGI and SG2 spectra. Since the positional stability 
was better when the SG2 spectrum was taken than the SGI spectrum for the observations 
of both positions, only the NG data taken simultaneously with the SG2 are used for the 
data at both positions. Then, the SGI spectrum is scaled to match with the SG2 spectra 
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in the spectral region of 7.3 to 7.9 /im. The scahng factors are 0.9 and 0.8 for Position 5 
and Position 8-1, respectively. The scaling of the SGI spectrum is taken into account in 
the derivation of the band intensity ratios (§3) and does not make serious effects on the 
following discussion. 



Results 



The resultant spectra toward our target positions are shown in Figure O The UIR 
bands are clearly seen at 3.3, 6.2, 7.7, 8.6 and 11.3 /i m in every spectruni . A weak feature 



20091 ). The hydrogen 



around 5.70 /xm is also seen at Positions 7 and 8-A (iBoersma et al. 
recombination line of Bra 4.05 /im is detected in four of the spectra, and the hydrogen 
recombination lines of Br/3 2.63 yum, Pf7 3.75 /im, and Pf/3 4.65 /im as well as the forbidden 
lines of [Ar II] 6.98 /im, [Ar III] 8.99 /im, [S IV] 10.51 /im and [Ne II] 12.81 /im are detected 
in three of the spectra. We note that a small bump seen around at 9.6 /im is an artifact, 
originating from the latent of the S9W exposure frame taken just before the SG2 exposure 
frames. Therefore, the spectral data from 9.4 to 9.8 /xm are not used in the following model 
fit, and we cut 9.4-9.8 /im from spectra in Figure [3l 



To derive the intensity of each UIR band and emission line, we fit the observed spectra 



with 



14 



11 



(A-AO^ + (7../2)^ 



5^ Cfc^ • exp 



1 (A-A,J^ 
2 '71/(8 -1112) J' 



fc=0 fc;=l 

where A is the wavelength. The first term represents the continuum, which is modeled with 
a polynomial function of the 5*^ order and is constrained to be non-negative. The second 
and third terms correspond to the UIR bands and the emission lines, respectively. For the 
UIR bands, we include 17 components centered at 3.30, 3.41, 3.46, 3.51, 3.5 6, 5.70, 6.22, 



6.69, 7.60, 7.85, 8.33, 8.61, 10.68, 11.23, 11.33, 11.99, 12.62 /im according to 



Smith et al. 
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( 120071 ) ■ Except for the 3.46, 3.51 and 3.56 /im components the UIR band components have 
band widths larger than or similar to the spectral resolution of AKARI/IKC and are thus 
modeled with Lorentzian profiles (the second term), where A^, is the center wavelength, 
7fcj is the FWHM and bk, is the height of each component. As for the 3.30, 3.41, 5.70, 
6.22, 6.69, 7.60, 7.85, 8.33, 8.61, 10.68, 11.23, 11.33, 11.99, 12.62 components, is fixed 
to the best-fit value obtained for the spectrum that has the highest S/N ratio with the 
spectral resolution of AKARI/IRC as the minimum value. The adopted value of 7^, for 
each component is summarized in Table [31 Only the height b^i is left as a free parameter. 
The integrated intensity of each component is calculated as nbki ' Ih/'^- following 
analysis, the 7.7 /xm band is defined as a combination of the 7.60 and 7.85 /im components, 
and the 11.3 /im band as a combination of the 11.23 and 11.33 /im components. We note 
that the 12.6 /im band is defined as one component of the 12.62 /xm co mponent, be c ause 



Smith et al. 



(120071), 



the red-wing component, such as the 12.69 yum component suggested by 
is not detected at a significant level in all the spectra due to the poor S/N and the low 
spectral resolution. 

The 3.46, 3.51 and 3.56 /im components as well as the emission lines, which have 
the band widths smaller than the spectral resolution of AKARI/IRC, are modeled with 
Gaussian profiles (the third term), where \kg is the center wavelength, -^^^ is the FWHM 
and Ckg is the height of each component. In the fit, 7^^ is fixed to match with the spectral 
resolution of AKARI/\YiC at the corresponding segment. Only the height Ckg is a free 
parameter in the fitting. The integrated intensity of each Gaussian component is given by 
(7r/ln2)-'^/^ • Ckg • 7kg/2. The adopted values of 7^^ for the 3.46, 3.51 and 3.56 jim components 
and the emission lines are summarized in Tables |4] and [5l respectively. 

The best-fit model spectra given by Eq. ([T]) are plotted together with the observed 
spectra in Figure HI The residual spectra are also plotted in the lower panel of each plot. 
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The derived intensities of the major UIR bands at 3.3, 6.2, 7.7, 8.6, 11.3, and 12.6/im 
and the emission hnes of Bra , Br/3 , Pf7 , Pf/3, [Ar II]6.98, [Ar III]8.99, [S IV]10.51, and 
[Ne II]12.81 are summarized in Tables E] and [TJ respectively, where those detected with more 
than 2(T are indicated. The uncertainties in the intensity are estimated from the fitting 
errors taking account of the observational uncertainties. 

As shown in Table [71 the 3.3, 6.2, 7.7, 8.6, 11.3, and 12.6 yum bands are detected at 
every target position except for the 12.6 fim band at Position 5, which shows the faintest 
emission. The 5.70 fim band are detected only at Positions 7 and 8-A. On the other hand, 
as shown in Table [71 all of the hydrogen recombination lines of Bra 4.05 yum, Br/3 2.63 fim, 
Pi'-f 3.75 yum, and Pf/3 4.65 fim and the fine structure lines from the ionized gas of [Ar II] 
6.98 fim, [Ar III] 8.99 fxm, [S IV] 10.51 fxm and [Ne II] 12.81 /im are detected at Positions 7, 
8-1 and 8-2, whereas none of them are detected at Positions 1, 2, 3, 4, 5 and 6 except for 
Bra detected barely at Position 5. Taking account of th e high ioniz ation potentials of 27.63 



19731 ) ■ respectively. Positions 



eV, 34.83 eV and 21.56 eV to form Ar^+, 8^+ and Ne+ flAUenl 
7, 8-1 and 8-2 are exposed to the hard incident radiation field powered by young massive 
stars and are associated with H II regions. This view is consistent with the radiation field 
conditions suggested by the IRAS and AKARI colors. Based on the characteristics of 
the observed spectra, we classify the targets into two groups: "Group A", which includes 
Positions 1, 2, 3, 4, 5 and 6 and "Group B", to which Positions 7, 8-1 and 8-2 belong. The 
members of Group A are supposed to be exposed to incident radiation fields of weak to 
moderate intensities and consist mostly of molecular clouds and PDRs. Group B members 
are all associated with H II regions. 

We investigate th e effects of extinction on the spectra observed at the present target 



positions according to 



Dobashi et al. 



(l2008[ l. The visual extinction Ay toward the present 



target positions ranges from 0.0 to 2.5 as shown in the last row of Table [T] We assume 
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the " LMC avg" extinction curve provided by IWeingartner fc Draind (120011 ) to estimate 
the infrared extinction. We also estimate the value of Ay from the observed ratio of Br/3 
to Bra at Posi tions 7, 8-1, and 8-2, ass uming the Case B condition of = lO^K and 



provided by 



lO^cm ^ (IStorev fc Hummer 



Dobashi et al. 



19951 ). The UIR band intensities corrected with A\ 



20081 ) differ from those corrected with Ay estimated from the 



Case B condition by l ess than 10 % at Posi tions 7, 8-1, and 8-2. We adopt those corrected 



with Ay provided by 



Dobashi et al. 



(j2008[ ) at these positions for consistency with the 
other target positions. The effect of extinction correction on the UIR band ratios is small 
(< 15%), and does not affect the following discussion. 

Next, we evaluate a contribution from the unresolved emission line PfS at 3.30 fim 
to the 3.3 fim band, and that from Pfa at 7.46 /im to the 7.7/im band. We assume that 
the intensities of PfS and Pfa are equal to 9.3% and 30.1% of that of Br/3, respectively, 
according to the Case B condition and subtract them from the intensity of the 3.3 fim and 
7.7 fim bands. In the present target positions, the contribution from Pf6 to the 3.3 fim band 
is less than 10% and that from Pfa to the 7.7 fim band is less than 3%, both of which are 
similar to the measurement uncertainties and thus do not affect the results. The intensity 
for which these corrections are applied is also listed in the lower row for each position in 
Tables O and [71 The values of the corrected UIR band ratios are listed in Table [SI The 
effect of these corrections on the values of the UIR band ratios is less than ~ 10% for all 
the targets. We estimate the uncertainties in /6.2/^m and /r.T^m as the difference between 
the intensities with and without scaling, which dominates over the fitting errors. 

Figures [5K-f show the plots of the corrected UIR band ratios of the /s.a^m/Ai.s/im, 

-^6.2^tm/Al.3Atm! h.7 nm/ hl.S fim, I 8. 6 fim/ hi .3 nm, h2. 6 fim/ 1 11. 3 fim, and /6.2 ^m/-^7.7 ^tm agaiust the 

IRAS color of h^fim/Iufim- Figures [5|;-1 plot those ratios against the AKARI color of 
-^L24/-^sii- Group A forms a sequence with a positive slope in the plots of h.s^m/ hi.Sfim, 
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h.2f,m/hi:sfim^ h.rf^m/hi.Sf^m and /g.eMm/Ai.SMm against the IRAS and ^iC^i?/ colors, while 
Group B does not follow the sequence in the plots of these band ratios. The /i2.6Atm/Ai.3/im 



and Jfi 



2^J.m/ h.7 



ratio does not show a systematic trend with the colors either for Groups 



A or B. Little variation fo und in the 6.2 um to the 7.7 fim band ratio is similar to the trend 



seen in external galaxies (IGalliano et al. 



20081 ) ■ The present observation shows that the 



3.3 /im band is weak in Group B compared to Group A. The results are discussed in the 
next section. 



4. Discussion 



4.1. Variation in the UIR band ratios 



Studies based on quantum chemical calculations as well as laboratory experiments show 
that the ionization fraction and the molecular size of PAHs are the prirnary factors to affect 



the relative intensi ty of the UIR band features (e.g. 



AUamandola et al. 



1999 



Draine fc Li 



2007 



Tielens 



20081). The ionization fraction of PAHs is determined by the balance between 
the photo-ionization and the recombination with ambient electrons and thus it is expected 

1 /2 

to be controlled by Go ■ , where Go is the intens ity of the radiation field in units of 



the solar vicinity value (1.6 x 10 ^Wm 



is the gas temperature (IBakes fc Tielens 



Habing 



1994 



larger fractions of positively ionized PAHs flBakes et al. 



1968h. is 



Bakes et al 



;he electron density, and Tg 



200 ll ). Higher Go/rif, ratios favor 



200ll ). Laboratory experiments and 



theoretical studies show that the ionization of PAHs enhances the intensity of the band 
featu r es in the 6-9 /xm region re l ative to the features in the 11-14 /xm region ( iDeFrees et al. 



1993 



Szczepanski fc Vala 



1993 



AUamandola et al. 



1999 



Bakes et al. 



20011 ). Therefore, 



the ionized-to-neutral band ratios (e.g. , Ie.2^^m/ hi.si^m, I^.^^,m/hl.3^lm, h.Q^lm/ hi.z^im) 
are supposed to indicate the ionization conditions of PAHs. Observational studies report 
variations in Ij.r^m/hi.Sfim and /s.e^jm/Ai.s^tm within a reflection nebula along the distance 
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stars with different spectral types fISloan et al. 



from the central star (jBregman fc Temil l2005l: jJoblin et al.lll996[ ). among Herbig Ae/Be 



of the star- for ming galaxy NGC 6949 (jSakon et al. 



20051). bet ween the interarm and arm regions 



outer Galaxy (jSakon et al. 



20071 ). and between the inner and the 



20041 ). They are reasonably interpreted by the difference in the 



ionization conditions of the band carriers. 

While very small PAHs (tt-c <~ 10^, where uq is the number of carbon atoms in 
the PAH) radiate strongly at 3.3 /im, large PAHs radiate mostly at longer wavelengths: 
PAHs as large as uq ~ 10^-10^ efficiently convert the absorbed energy t o the 6.2, 7.7 
and 8.6 /im bands and P AHs as large as uq ~ 4000 to the 11.3 /im band (jSchutte et al. 



1993 



Draine fc Li 



20071 ). The 6.2 /im and 7.7/im bands are attributed to stretching modes 
of C-C bonds, while the 3.3/im, the 8.6 /im and the 11.3 /im bands are attributed to 
stretching m odes, in-plane bending m odes, and out-of-plane bending modes of C-H bonds, 



respectively (lAUamandola et al 



19891 ). Therefore the ratios of the short-to-long wavelength 



UIR bands from the same bonds (e.g., Is.s^m/h 



/^m/ 11.3 /xni) -'6.2 /xm/-^7 



the average temperature or the size distribu tion of PAHs (e.g. 



7,,,rn) can be used to infer 



Jourdain de Muizon et al. 



1990 



Sales et al. 



2010 



Boersma et al. 



2010l ). The NIR to MIR spectra we discuss here are 



obtained from the same region in the sky and thus we can discuss the band intensity ratios 
of the major UIR bands at 3.3, 6.2, 7.7, 8.6 and 11.3/xm concurrently. 

As described in |J3|, Group A forms a sequence on the diagrams of the IRAS and 



AKARI colors v.s. the ionized-to-neutral UIR band ratios: /6.2^m/Ai 



^m/ 11.3 /im; 7.7 /im/-^11.3 /jm 

and /8.6Atm/Ai.3^m- The sequence suggests that a larger fraction of PAHs is ionized as 
the radiation field becomes stronger. However, Group B, whose radiation fields are much 
stronger and harder than those of Group A, does not follow the sequence. This can be 
attributed to the lower ionization fraction of PAHs due to an increase in the recombination 
under the high electron density in H II regions relative to molecular clouds or PDRs. 
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Papoularl ( 120051 ) argues that hydrogen impact might also play a role in the excitation of 
PAHs in PDRs and molecular clouds. 

There is a positive correlation between IRAS and AKARI colors and the short-to-long 
wavelength UIR band ratio, /s.a^m/Ai.SAtm in Group A. This can be interpreted in terms 
of an increase in the excitation temperature of PAHs with the IRAS and AKARI colors 
in Group A. Destruction of PAHs is expected to be inefficient in environments of Group 



A (iMicelotta et al 



2010a b 



20111 ) and the size distribution of PAHS does not change 
considerably. Thus the excitation temperature of PAHs is mainly controlled by the hardness 
of the incident radiation field, but not by its intensity. On the other hand, the IRAS 



and colors indicate its intensity, but not the hardness directly (jSakon et al. 



2006 



Onaka et al. 



2007bl ). The correlation seen in Figure [5] thus indicates that the incident 
radiation field becomes harder as the intensity becomes larger in Group A. This is a 
reasonable consequence of strong incident radiation fields, for which a contribution from 
young massive stars becomes dominant. The present observations indicate this trend 
explicitly based on the NIR to MIR UIR band ratios. Harder incident radiation fields 
also increase the /6.2/.m/Ai.3Mm, h.r f,m/ hi/i/im, and hsiMm/ hi.si^m- We show in §4.2. 2| that 
the ionization fraction is the major factor for the increase in these ratios with a minor 
contribution from the hardness of the incident radiation field. 

Group B does not follow the trend of Group A and shows the weaker 3.3 yum band. It 



can be attributed to destruction of very small PAHs (nc 10^) inside 



are supposed to dominantly co ntribute to the 3.3 /xm band (ISchutte et al 



20071). 



Micelotta et al. 



i II regions, which 



1993 



Draine&Li 



(j2010af ) suggest that electron collisions dominate in the destruction 
of PAHs in a hot plasma for gas temperatures between ~ 3 x 10"^ and ~ 10^ K, where smaller 
PAHs are more efficiently destroyed than larger ones. The threshold temperatures depend 
sensitively on the assumed electronic dissociation energy. The difference of /s.s/im/Ai.s^tm 



- 19 - 



between Group A and B may thus result from the destruction by electron collisions, if it 
is a dominant process of PAH destruction in Group B. The absence of a systematic trend 
with the IR colors may be partly due to the small number of the data in Group B. 

4.2. Diagnostic diagram of the radiation field condition 

4-2.1. The diagnostic diagram of the UIR band ratios 

According to the discussion in the previous subsections and the results presented in ^ 
we explore possible diagnostic diagrams to investigate the physical conditions of the ISM 
by means of the UIR band ratios. The ratio /3.3^m/-^ii.3Atm is found to be a good indicator 
of the size distribution of PAHs, whereas /6.2Mm//ii.3A»m, h.i tim/ hi.?.^im, or h.&tixn/ hi.^^ira 
has been indicated to be a measure of the ionized fraction of PAHs as discussed above. 
The 8.6 /im band is typically weak and situated on the shoulder of the strong 7.7/im band 
due to the limited resolution provided by the IRC. Thus h.Q^ra/ is less reliable than 
the other two ratios and will not be considered. Taking account of these, we investigate 
two-band-ratio plots of h.-i^im/ hi.^i^ira v.s. h.r ^,m/ hi.3f,m or h.2^ira/ hi.^^ira as shown in 
Figures |6^ and b. 

In Figure |6t where /3.3^m/Ai.3/i,m is plotted against h.-j ^ra/ hi.ZiJ.ui, Group A forms a 
sequence from the bottom-left to the top- right. This sequence can be interpreted in terms 
of the change in the spectrum of the incident radiation field together with the change in the 
ionization fraction of PAHs. Group B, however, does not follow the sequence and is located 
distinctly from Group A. FigureEb shows a plot of /3.3^m/-^ii.3^fm v.s. /6.2^jm/-^ii. 3^*111- Group 
A forms a sequence similar to the plot of /3.3^m/Ai.3^jm v.s. /7.7;jm/Ai. 3^*111 and Group B 
data are also separated from Group A. The separation of Group B from Group A is more 
distinct in Figure than in Figure Eb- Lower left points of Group B in Figure may 



- 20 - 



overlap with the sequence of Group A if the sequence continues down to the lower ratios 
and would become difficult to be distinguished from Group A. The difference between the 
two plots comes from a slight decrease in the 6.2 /im band intensity relative to the 7.7/im 
band in Group B, in particular in Position 8-1 and 8-2 spectra. While the h.2^lm/ h.i ij.ni 
ratio of Group A does not show a clear dependency on the intensity and/or hardness of the 
radiation field (Figures and j), the average value of the lQ.2nm/ h.i iira ratio of Group B 
(0.303±0.053) is slightly smaller than that of Group A (0.375±0.051). This may suggest 



that even mediurn - sized PAHs are d estroyed in H II regions in Group B ( ISchutte et al. 



Draine fc Li 



2001 



Sales et al 



1993 



2010l ). Alternatively, de-hydrogenation might also explain 



this trend. Since C-H in-plane bending modes contribute to part of the 7.7/im band, H-loss 
reduces the intensity of the 7.7/im band. However, there is a large scatter in /6.2/im/-^7.7^tm 
among each group. We cannot rule out a possibility that there are other causes for the 
variation in lf,.2^xalh 

In the present results. Groups A and B can be distinguished more clearly in the 
diagram of /3.3A»m//ii.3/im v.s. Ir.rfim/hi.atim than that of h.^ii^nj hi.'i^m v.s. Ie.2tMm/hi.3tim- 
It suggests a potential of the diagram of Is.Sf^m/ hi.Sfim v.s. Ij.t ^m/ hi.znm as a diagnostic 
tool for the radiation field conditions. The 6.2 /im band is usually thought as a more reliable 
indicator than the 7.7/im band, since the 6. 2 /xm consists of a s i ngle comp orient, w hile the 



7.7/im contains more than one component ( iPeeters et al. 



2002 



Smith et al. 



20071 ). The 



present results indicate that the 7.7/im band might be a better indicator than the 6.2 /im 
band as far as the ionization fraction of PAHs are concerned probably because it is less 
sensitive to the variation in the size distribution. The number of Group B targets is small, 
however, and thus we need further investigations to ensure if this conclusion can be applied 
to the general ISM with a wide range of physical conditions. 



Vermeij et al. 



(I2OO2I ) investigate the variations in the MIR UIR bands among the 
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H II regions of the LMC based on ISOPHOT/PHT-S observations and find that the 
h.2^lu^/ hl.2,^ln'^ and Ij.t nm/ hi.Sfim ratios are systematically smaller in HII regions in 30 
Doradus and the Small Magellanic Cloud (SMC) than in Galactic and non-30 Doradus 
regions. They suggest that the segregation may be attributed to the difference in the 
molecular structure of PAHs, proposing that compact PAHs may dominate in the SMC 
and 30 Doradus regions. The present targets do not include 30 Doradus itself. Also none 
of them have the incident radiation fields of the intensity similar to the 30 Doradus region. 
The 11.3 /im band is due to solo C-H out-of-plane bending modes and probes long straight 
edges, whereas the 12.6 /im band is due to trio C-H out-of-plane modes and probes corners. 
The contribution from compact PAHs enhances the 11.3/im band relative to the 12.6 /im 
band. The different compactness of PAHs is supposed to appear more distinctly in the 
h2.6fj.u1 /hi. 3 fim ratio than in the /s.s^m/Ai.s^tm- The ratio of h2.6fj.n1/ hi. 3 fim is 0.40 ±0.10 
and 0.44 ± 0.03 for Group A and B, respectively. Although the intensity of the 12.6yum 
band has a large uncertainty, the ratio does not show a distinct difference as seen in 
h.3fj,ni/hi.3tim, Suggesting that the molecular structure does not change appreciably among 
the present targets and is not the major factor for the difference in /s.s/xm/Ai.SAtm- 



Recently 



Hony fc SAGE-Spec Team 



(120111 ) also investigate the variations in the MIR 



UIR bands based on Spitzer/IKS observations of the LMC. Their targets include molecular 
clouds, PDRs and H II regions as in the present study, but the number of the targets is much 
larger. They find that the ratio h.7 /im/ hi.3fjm is mostly controlled by the ionization fraction 
of PAHs and weakly depends on the [Ne II] /[Ne III] line ratio, which they interpreted as 
the indicator of hardness of the incident radiation field. The present results show a similar 
trend, but also suggest a systematic trend with the hardness, once molecular clouds and 
PDRs are separated from H II regions by use of the ratio /a.s^m/Ai.a^im- The sequence seen 
for Group A suggests that the excitation of PAHs as indicated by h.3fj.m/hi.3fj.m is enhanced 
with the ionization fraction, which is indicated by h.r fjm/ hi. 3 fjm or /6.2/^m/Ai.3/^m- Thus the 
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ratio /3.3^m/Ai.3/im IS a useful measure for the incident radiation field conditions for Group 
A targets. This interpretation suggests that the data at the lower left in Figure E] are objects 
in an early stage of the cloud evolution and those at the upper right are more evolved 
PDR-type objects. Further observations are important to confirm the interpretation. In 
both plots of Figure [6l the ratio /a.s^m/Ai.s^tm plays a crucial role not only to separate 
Group A from B, but also to estimate the incident radiation field conditions of the target 
in Group A. 



4-2.2. Comparison with model 

The previous subsection suggests possible diagnostic diagrams based on the UIR band 
ratios and qualitative interpretation is given. In this subsection we employ simple models 
of PAH emission and investigate the observed diagrams quantitatively. 

To derive the intensity of the UIR bands emitted from a mixture of neutral and 
ionized PAHs of v arious sizes, we emp loy a simple theoretical model of infrared PAH 



emission following 



Schutte et al. 



( I1993I ). The A-coefficients of the UIR bands in 3-20 /im 



are calculat ed from the infrare d cross sections of neutral and ionized PAHs recently 



pr ovided by 



by 



Draine fc Lil (120071 ) and the in ternal energy to temperature relation provided 



Draine fc Lil (l200l[ ) is adopted (see also 



Boersma et al 



2010 



Bauschlicher et al 



2OIOI) 



The number of carbon atoms in PAHs is assumed to be distributed between n™™ and n™^^, 
and the size distribution of PAHs is assumed to be given by the same power-law function 
as graphite grains. Details on the model calculation are given in Appendix lAl 

The UIR band ratios of /s.s^m/Ai.s/im and h.r ^m/hi.s^im are calculated for different 
temperatures of the heating source and ionization fractions of PAHs /i. Two sets of 



'-c ' "-c 



are calculated in the following analysis: Case I assumes (n™™, n™' 



(20, 
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4000). The UIR band ratios of /s.s^tm/Ai.a^tm and Ij.r^m/hi.s^im calculated for various 
and /i are plotted in Figure [7^. Compared with the locations of the Group A data points 
on the plot, the bottom-left edge of the sequence corresponds to T^, ~7500K and /i ~30%, 
while the top-right end of the sequence to that with ~12000K and /i ~60 %. The 
range of the effective temperature indicated is compatible with the effective temperature 
of A-type main sequence stars of 7000-lOOOOK. Thus the model calculations support the 
interpretation in §4.11 and the observed sequence of Group A can be accounted for by the 
increase of the effective temperature of the heating source, which increases the excitation, 
accompanied by the increase in the ionization fraction of PAHs without the change in the 
size distribution. The observed data of Group B, on the other hand, are distributed in 
a region corresponding to T* ~6500-8000 K and /i ~40-60 % for Case I, which is totally 
inconsistent with the fact that the radiation field of Group B should be much harder than 
that of Group A. 

Case II assumes (ng''^, n^^^) = (100, 4000). The UIR band ratios of h.3^,m/hl.3^^m 
and Ij.T ^m/ hi.3fj.m calculated for various T* and /i are plotted in Figure [7b. The observed 
data of Group B are distributed in a region corresponding to T^, ~30, 000-40, 000 K and 
/i ~20-40 %, which is compatible with the effective temperature of 0-type main sequence 
stars of 30,000-50,000 K. The small /i is consistent with the increase in recombination with 
electrons in high electron density environments of H II regions. Therefore, the deviation of 
Group B from the sequence of Group A on the plot of /s.s^m/Ai.SAtm v.s. Ij.j hi.Sfim is 
reasonably accounted for by the difference in the minimum size of PAHs. 

Therefore, the model calculations quantitatively support the interpretation of the 
trend seen in the plot of /3.3^m/-^ii.3Atm v.s. h.r ^m/ hi.Sfim- We stress here that an accurate 
measurement of I3.3 fim/ hi.s fim ratio should be the key to obtain the average excitation 
temperature of PAHs, which is determined by the size distribution of PAHs and the 
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hardness of the incident radiation field. The present calculation of Case I and II is made 
only for representative purposes and the choice of n™™ is rather arbitrary. The value of 
-^s.SAtm/Ai.s/im IS scusitive to the specific heat as well as the size distribution of PAHs that 
we assume in the calculation. The actual values of T*, and n™"^ inferred from the plot 
depend on the assumed PAH properties such as the size distribution and the specific heat, 
which are also questionable. 



5. Summary 

We present the results of NIR to MIR slit spectroscopic observations of the diffuse 
radiation toward nine positions with different radiation field conditions in the LMC with 
AKARI/IKC. We obtain continuous spectra from 2.55 to 13.4yum of the same slit area, 
which allow us to investigate variations in the relative intensity of the UIR bands from 2.55 
to 13.4 /im emitted from exactly the same region. 

The target positions are selected based on the IRAS colors of I25 ^im/Iufira and 
hofim/hoofimj which indicate star formation activities, to cover a wide range of the intensity 
of the incident radiation field. The AKARI coloi of /l24/-^sii at the present target positions 
shows a similar trend to that of the IRAS l25fj.m/ Ii2fj.m color, confirming that the selection 
based on the IRAS colors is in fact relevant to the purpose of the present study. 

A series of the major UIR bands at 3.3, 6.2, 7.7, 8.6 and 11.3 fim are clearly detected 
in every obtained spectrum, and the hydrogen recombination lines of Bra 4.05 /im, Br/3 
2.63 /im, Pf7 3.75 /im, and Pf/3 4.65 fim and the fine structure lines of [Ar II] 6.98 /im, 
[Ar III] 8.99 fim, [S IV] 10.51 fim and [Ne II] 12.81 fim, signatures of the presence of ionized 
gas, are detected in some of them. According to the ionization gas signatures in the 
spectra, we classify the present nine target positions into two groups: those without the 
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strong ionized gas signatures (Group A) and those with the signatures (Group B). Group A 
positions are supposed to be in relatively quiescent radiation field environments and those 
of Group B are in harsh radiation field environments powered by young massive stars. This 
view is consistent with the radiation field conditions suggested by the IRAS and AKARI 
MIR colors. 

Group A shows a sequence on the plots of the UIR band ratios of /s.s^m/Ai.a^im, 
/6.2Mm//ii.3/i,m, h.j fim/ hi.3 and h.s fim/ hi.3 agaiust the IRAS and AKARI colors, but 
Group B does not follow the sequence. These results can be interpreted in terms of the 
facts that (1) in Group A, PAHs are heated to higher excitation temperatures and their 
ionization fraction increases as the radiation field becomes harder and stronger and that (2) 
in Group B, very small PAHs (nc < 100) are efficiently destroyed, possibly due to electron 
collisions, and the ionization of PAHs is suppressed by an increase in the electron density 
inside H II regions. The present observations also show that the incident radiation field 
becomes harder as the intensity increases in Group A based on the UIR band ratios and the 
IR colors. There is little variation in I6.2fim/l7.7 fim as reported in previous studies and we 
find no systematic trend against the colors. 

The observed data points of Groups A and B are well separated on the plot of 

I3.3fj.m/In.3fim V.S. /y.? ^m/ Al.3 aS Well aS that of J3.3^m//ll.3Mm V-S. h.2 fim/ hl.Z fim- TheSe 

trends can be interpreted in the same way as described above, suggesting a potential of the 
diagram of /a.s^m/Ai.s^tm v.s. /T.r^m/Ai.s^tm as a useful diagnostic tool for the radiation field 
conditions. Simple model calculations support the interpretation quantitatively. Further 
investigation is needed to ensure the applicability to this diagram a wide range of objects. 

The present study shows the importance of the UIR band at 3.3 /im. The ratio 
h.z lira/ 1 11. 3 lira plays a crucial role not only to separate Group A from B, b ut also to es t imate 



the incident radiation field conditions of the target in Group A. Recently 



Seoketal 



mm 
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report the detection of the 3.3 fim UIR band at the supernova remnant (SNR) N49 in 
the LMC, suggesting the presence of significant processing of PAHs in SNR shocks. The 
3.3 fim band provides significant information on the size distribution and/or the excitation 
conditions of PAHs. 

This work is based on observations with AKARI, a JAXA project with the participation 
of ESA. The authors thank all the members of the AKARI project and the members of the 
Interstellar and Nearby Galaxy team for their help and continuous encouragements. The 
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AKARI. We thank A. Kawamura and Y. Fukui for providing us the LMC 12C0 survey 
data at 2.7 mm taken by the NANTEN millimeter- wave telescope of Nagoya University. 
This work is supported in part by a Grant-in- Aid for Scientific Research from the Japan 
Society of Promotion of Science (JSPS). 



A. Appendix 



The present model calcula tion basically 



"ollows 



model parameters provided by 



Boersma et al. 



2010 



Bauschlicher et al. 



Draine fc Lil (1200 ih and 



Schutte et al. 



1993) with the recent 



Draine fc Lil (120071 ) (see also 



20101). The emission intensity due to an IR active 



fundamental vibrational transition, i, from level (v-l) to level v, in a j-type PAH molecule 
with a total internal vibrational energy E is given by 



I{j, E, i, v) = huivA 



(Al) 



where h is the Plank constant, v is the vibrational quantum number, i/j is the frequency 



of the emitting mode, is the Einstein coefficient of the 1^0 transition, Pj{E) is the 



-27- 



total density of vibrational states at total energy E, i.e., the number of ways the energy E 
can be distributed over all available states, and Pj,r{E — vhvi) is the density of vibrational 
states for all modes except the emitting mode at a vibrational energy E — vhvi. For the 



Einstein coefficient of the 1 — ?■ transition, A^'^ is given by 



A 



1,0 
hi 



Svrc 

^4 ^j,int,i) 



(A2) 



where Aj is the wavelength of the emitting mode, CTj,int,i is the cross section of the i-th mod e 



i ntegr ated over the wavelength. We adopt the values described in Table 1 of 



Draine fc Li 



(120071 ) as (Tj,int,i for each type of PAH molecules. Here we calculate the model spectra with 
various parameters to semi-quantitatively compare with the observations and thus adopt a 
simple thermal approximation. The validity of the thermal approximation has been studied 
to a large extent and it is shown that the presence of the size distribution alleviates the 



difference and the effect on the relative band intensities is smal 



study (lAUamandola et al. 



1989 



Schutte et al. 



1993 



Draine fc Li 



enough for the present 



20011). 



In the thermal approximation, the emitted intensity of a j-type molecule with internal 
energy E in the i-th mode, from level v to level (f-1), is described by 



I{j,E,i,v) = huivA-'- exp{-vhi'i/ kTj{E))[l - exp{-hi'i/kTj{E))] \ 



(A3) 



where Tj{E) is the vibrational excitation temperature of a j-type molecule with E and 
where k is the Boltzmann constant. The sum of equation (IA3p from f = 1 to f = oo, i.e., 
the total emitted intensity in the i-th mode, is given by 



Iij,E,t) = hu,A]l^[exp{hu,/kT,{E)) - 



(A4) 



In this approximation, the energy-temperature relation for a j-type molecule is given 



by 



(A5) 



i=l 
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where s is the number of vibrational modes of a j-type molecule equal to 3 ti mes of atoms 
minus 6, i.e., Sriatom ~6. We derive this relationship, using equations (2)-(8) in 



Draine&Li 



( 1200 ll ) for each type of PAH molecules. 



The total energy emitted in the i-th mode following the absorption of a UV/visual 
photon of frequency of u, f{j, u, i) is given by 

where J(j, i) / E) is the fraction of the total IR intensity emitted by a j-type molecule 
with the internal energy E in the i-th mode. Hence, the total emitted intensity in the 
z-th mode of a j-type molecule exposed by a star, whose spectrum is approximated with a 
blackbody, P(j, T*,«) is calculated as 

T., z) = r ""^'^f^^^h u, V, i)dv, (A7) 
where aj^u is the UV/visual absorption cross s ection of ?"-type mole cule, as which we adopt 



the values described in equations (17)-(20) in 



Schutte et al. 



fll993[ ). Then, the flux in the 



z-th mode from a certain size distribution of interstellar PAHs, exposed to a star with T*, is 
given by 

F(T„z) = ^npAH(j)^(j,T„2), (A8) 

3 

where tt-pahO) is the number density of the j-type PAH molecule. 

In order to evaluate the effect of depletion of very small PAHs quantitively, we calculate 
the ratio of F(T^,,i) relative to the FiT^.V), the model ratio of the emitted intensity in the 
z-th mode relative to that in the z'-th mode in two cases; Case I and Case II (see text). In 
both cases, we assume that the number density of interstellar PAHs is given by a power-law 
distribution with the size, i.e., the number of carbon atoms contained by PAH molecules, 
uq. If the number of PAH molecules per interstellar H atoms with a radius a between 
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a + da is proportional to a " and uq (x a* , the number of PAH molecules per interstellar H 
atoms with a number of carbon atoms uq between uq + dnc is given by 

dNpAH 



1 — a — 7 



N. 



(A9) 



H 



We adopt q;=3.5 and 7=8 as assumed in 
given by 



Schutte et al. 



(Il993l ). Then the number density is 



(AlO) 



AUamandola et al. 



(1l989r ) suggest 



In Case I, the minimum size of PAHs is set to nc=20. 
that smaller PAHs with ric < 20 are photolytically unstable. In Case II, we set the minimum 
size of PAHs nr=1 00. In both cases, the maximum sizes of PAH molecules is fixed at 4000. 
Draine fc Lil (120071 ) suggest that larger PAHs with nc > 4000 do not contribute to the UIR 
band features at 3-11 yum. Then, w e run the model t\Y ice with the cross sections for neutral 



PAHs and ionic PAHs provided by 



Draine Sz Lil (120071 ). respectively, and calculate the ratio 



of F(T^,,i) from a summation of the spectra of both components by varying the ionization 
fraction. 
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Table 1. Observation log and parameters 



Position ID 


Date 


Obs.ID 


AOT 


Disperser 


Center position 


of the slit 














«2000 


<520()0 




Position 1 


2006 October 14 


1400318.1 


IRC04 b:Ns 


NG, SGI and SG2 


05'"38™20.=4 


-70°07'25" 


1.89 


Position 2 


2006 October 19 


1400330.1 


IRC04 b:Ns 


NG, SGI and SG2 


05'^ 39™ 43. =7 


-70°00'34" 


0.721 


Position 3 


2006 October 2 


1400346.1 


IRC04 b:Ns 


NG, SGI and SG2 


05'M8™08.=2 


-69°52'52" 


1.90 


Position 4 


2007 May 15 


1402426.1 


IRC04 b:Ns 


NG, SGI and SG2 


05''25™55.=4 


-66°10'30" 


0.257 


Position 5 


2006 October 16 


1400324.1 


IRC04 b:Ns 


NG, SGI and SG2 


0511430142 sg 


-69°22'32" 


0.826 


Position 6 


2007 May 18 


1402422.1 


IRC04 b:Ns 


NG, SGI and SG2 


05'^26™06.=3 


-68°36'04" 


0.630 


Position 7 


2006 October 16 


1400334.1 


IRC04 b:Ns 


NG, SGI and SG2 


05'^39™17."67 


-69° 30' 14" 


0.541 


Position 8 


2006 October 16 


1400320.1 


IRC04 b:Ns 


NG, SGI and SG2 


05'^39™57.=5 


69°45'27.5" 


2.11 


LMC-off 


2006 November 10 


1500719.1 


IRC04 a:Ns 


NP, SGI and SG2 


06'"00™00.=0 


-66°36'30" 




LMC-off 


2006 November 10 


1500720.1 


IRC04 a:Ns 


NP, SGI and SG2 


06'" 00™ 00. =0 


-66°36'30" 





^From 



Dobashi et al 



1 20081) 



Table 2. IRAS and AKARI colors^ 



Position ID 


l25pt.ni/ -fl2fim 


-'^eOfimZ-flOO/jm 


I'L2i/ -fsil 


Position 1 


l.OlltO.Ol 


0.371±0.004 


0.58±0.04 


Position 2 


1.17±0.01 


0.359±0.001 


0.76± 0.05 


Position 3 


1.30±0.05 


0.317±0.004 


1.3±0.1 


Position 4 


1.47±0.05 


0.435±0.004 


1.4±0.1 


Position 5 


1.67±0.04 


0.474±0.006 


1.5±0.1 


Position 6 


1.82±0.02 


0.467±0.004 


l.SibO.l 


Position 7 


4.43±0.16 


0.751±0.014 


6.9±0.5 


Position 8-1 


4.44±0.34 


0.730±0.010 


7.9±0.5 


Position 8-2 


4.44±0.34 


0.730±0.010 


6.9±0.5 



''The values are the ratios of the radiances in W m ^ Hz ^ 
arcsec"^ . 
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Table 3. UIR band parameters modeled by Lorentzian profiles 





7fei [/^m] 




n r\AO 

U.U4Z 




U.UDU 


0. / u 


u.zu 


6.22 


0.20 


6.69 


0.48 


7.60=^ 


0.34 


7.85=^ 


0.44 


8.33 


0.50 


8.61 


0.34 


10.68 


0.22 


11.23"^ 


0.21 


11. 33*= 


0.38 


11.99 


0.54 


12.62 


0.53 



^the components of the 
7.7 /im complex feature. 

'^thc components of the 
11.3 /im complex feature. 



Table 4. UIR band parameters modeled by Gaussian profiles 



^kg [t^ra] 7fc^ [^J.ra] 



3.46 
3.51 
3.56 



0.034 
0.034 
0.034 
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Table 5. Emission line parameters modeled by Gaussian profiles 



Line 




7feg [/^m] 


Br/3 


2.63 


0.034 


Bra 


4.05 


0.034 


Pf7 


3.74 


0.034 


Pf/3 


4.65 


0.034 


[Ar II] 


6.98 


0.114 


[Ar III] 


8.99 


0.200 


[S IV] 


10.51 


0.200 


[Nc II] 


12.81 


0.200 
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Table 6. Observed and corrected intensity of the UIR bands^ 



Position ID 


^3 . 3 (.mi 


^6.2(im 






.fll.SMm 


-fl2.6nm 


JT UoililOli ± 


n fififl-l-n 049 


9 1 9+n 1 Q 






'\ 44-1-0 9K 

O.^'i^U. ZO 


1 1 1 +n '^Q 




U. t OO^U.UIU 








Q 074-0 ^9 


1 9n+n 4^ 


Position 2 










9 SO-t-0 94 


1 OA^n '^^ 
1 .u^^u. ou 




0.531±0.065 


1.98±0.18 


4.25±0.28 


0.896±0.130 


2.96±0.26 


1.08±0.36 


Position 3 


0.980±0.041 


3.40±0.22 


9.92±0.34 


1.98±0.16 


4.89±0.32 


2.95±0.45 




1.09 ±0.04 


3.56±0.23 


10.5±0.3 


2.32±0.18 


5.66±0.37 


3.21±0.49 


Position 4 


0.693±0.048 


2.87±0.15 


7.15±0.23 


1.75±0.10 


3.27±0.21 


1.41±0.30 




0.704±0.048 


2.89±0.15 


7.21±0.24 


1.79±0.10 


3.33±0.22 


1.43±0.31 


Position 5 


0.357±0.043 


1.32±0.14 


3.81±0.25 


1.01±0.11 


1.68±0.24 






0.370±0.049 


1.35±0.14 


3.89±0.26 


1.08±0.12 


1.79±0.26 




Position 6 


0.858±0.035 


3.11±0.15 


8.36±0.25 


1.79±0.11 


2.99±0.22 


0.99±0.32 




0.891±0.036 


3.15±0.15 


8.53±0.25 


1.88±0.11 


3.13±0.23 


1.02±0.33 


Position 7 


4.85 ±0.04 


23.1±0.2 


63.7±0.3 


13.0±0.1 


25.3±0.3 


11.92±0.45 




4.81 ±0.05 


23.4±0.2 


64.2±0.3 


13.6±0.1 


26.4±0.3 


12.21±0.46 


Position 8-1 


4.07 ±0.07 


14.5±0.2 


52.3±0.5 


13.4±0.2 


28.3±0.5 


15.04±0.74 




4.44 ±0.09 


15.3±0.2 


55.4±0.5 


16.0±0.2 


33.2±0.6 


16.52±0.82 


Position 8-2 


1.13 ±0.10 


4.76±0.39 


17.8±0.6 


4.79±0.27 


7.01±0.57 


3.05±0.83 




1.16 ±0.11 


5.01±0.41 


18.6±0.6 


5.70±0.33 


8.24±0.67 


3.35±0.91 



Note. — The UIR bands detected with more than 2 a are Usted. For each position, the upper row 
indicates the original band intensity and the lower row shows the intensity for which the extinction 
and the contribution from the hydrogen recombination lines are corrected (see text). 

^In units of 10~^*W m~^ arcsoc"^. 



- 39 - 



Table 7. Observed and corrected intensity of the emission lines' 



Position ID -fBr,82.63 -fBrQ4.05 ^Pf73.74 -fpf^4.65 -'^[ArII]6.98 ^[ArIII18.99 ^[SIVllO.51 ^[NcII]12.Sl 

Position 1- - - - - -- - 

Position 2 - - - - ____ 

Position 3 - - - - - _ 

Position 4 - - - - ____ 

Position 5 - 0.0599±0.0159 - - . . . _ 

0.0610±0.0160 - - - - - _ 

Position 6 - - - - ____ 



Position 7 


1.02±0.03 


2.04±0.01 


0.261±0.017 


0.457±0.022 


1.65±0.07 


1.48±0.05 


1.56±0.05 


3.94 ±0.12 




1.07±0.03 


2.08±0.01 


0.267±0.018 


0.464±0.023 


1.67±0.08 


1.57±0.05 


1.65±0.05 


4.03 ±0.12 


Position 8-1 


0.832±0.053 


1.70±0.02 


0.285±0.034 


0.279±0.029 


0.780±0.105 


1.99±0.08 


0.948±0.082 


3.63±0.20 




0.996±0.064 


1.85±0.03 


0.315±0.037 


0.297±0.030 


0.812±0.109 


2.51±0.11 


1.18±0.10 


3.97±0.22 


Position 8-2 


0.631±0.067 


1.27±0.03 


0.141±0.036 


0.262±0.042 


0.482±0.145 


1.39±0.09 


1.26±0.09 


1.83 ± 0.23 




0.755±0.080 


1.38±0.03 


0.156±0.040 


0.279±0.045 


0.503±0.151 


1.75±0.12 


1.57±0.11 


2.00 ± 0.25 



Note. — Emission lines detected with more than 2 a are listed. For each position, the upper row indicates the original band intensity 
and the lower row shows the extinction-corrected intensity (see text). 

''In units of 10~"'^*W arcsec"^. 
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Table 8. The UIR band intensity ratios 



Position ID 


-^3 . 3 fim / -^11.3 i-im 


l6.2 jjm/-fll.3 Mm 


h.T nm/Ill.'A Mm 


Is.e iJ.in/Ill.3 Mm 


-'^12.6Mm/-fll.3 Mm 


l6.2 Mm/ It.T i.in\ 


Position 1 


0.185±0.019 


0.559±0.069 


1.68±0.16 


0.409±0.052 


0.304±0.111 


0.331±0.035 


Position 2 


0.179±0.027 


0.671±0.085 


1.43±0.16 


0.303±0.051 


0.365±0.129 


0.466±0.053 


Position 3 


0.193±0.015 


0.629±0.058 


1.86±0.13 


0.411±0.042 


0.567±0.095 


0.337±0.025 


Position 4 


0.211±0.020 


0.867±0.073 


2.16±0.15 


0.537±0.048 


0.428±0.097 


0.401±0.025 


Position 5 


0.205±0.040 


0.752±0.136 


2.16±0.38 


0.601±0.112 




0.347±0.044 


Position 6 


0.284±0.024 


1.01 ±0.09 


2.71±0.21 


0.591±0.057 


0.326±0.108 


0.370±0.021 


Position 7 


0.182±0.003 


0.887±0.013 


2.42±0.03 


0.514±0.009 


0.461±0.018 


0.365±0.003 


Position 8-1 


0.133±0.003 


0.461±0.096 


1.66±0.33 


0.481±0.012 


0.496±0.026 


0.276±0.006 


Position 8-2 


0.141±0.018 


0.608±0.069 


2.26±0.20 


0.692±0.068 


0.406±0.115 


0.268±0.023 



Note. — The band intensity ratio for which the extinction and the contribution from the hydrogen recombination lines are 
corrected. 
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Fig. 1. — The sht positions of each pointed observation are overlain on the artificial 3-color 
image of the LMC, where the intensity of S7 (7/im) is indicated in blue, Sll (11 /im) in 
green, and L24 (24 /im) in red. The white crosses indicate the center positions of the slits. 
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(a) (b) 




R.A. (J2000) 

Fig. 2. — The enlarged slit position images of (a)Position 1, (b)Position 2, (c)Position 3, 
(d)Position 4, (e)Position 5, (f)Position 6, (g)Position 7, and (h)Position 8 are overlain on 
the Sll images of the LMC The white boxes indicate the slit of 1' length by 5" width in 
each observation. 
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R.A. (J2000) 



Fig. 2. — Continued. 
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CO 

PAH main | | | | 1 

sub HI I I I I LJI 

3 4 5 6 7 8 9 10 11 12 13 

Wavelength [|jm] 

Fig. 3. — Observed spectra of nine regions in the LMC; (a) Position 1, (b) Position 2, (c) 
Position 3, (d) Position 4, (e) Position 5, (f) Position 6 (g) Position 7, (h) Position 8-1, 
and (i) Position 8-2. The bottom panel indicates the positions of the UIR bands and the 
emission hnes. 
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Fig. 3. — Continued. 
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Fig. 4. — Modeled spectra are overlain on the observed spectra of nine regions in the LMC; 
(a) Position 1, (b) Position 2, (c) Position 3, (d) Position 4, (e) Position 5, (f) Position 6 
(g) Position 7, (h) Position 8-1, and (i) Position 8-2. The red solid lines show the best-fit 
model spectra (see text). The light-red dotted lines indicate the continuum component of a 
quintic function. The light-green, dark-green, yellow, coral, dark-pink, pink, light-blue and 
dark-blue lines show the UIR 3.3 /im, 3.4 /im, 6.2 /im, 7.6 /im, 7.8 /im, 8.6 /im, 11.2 /xm and 
11.3 /im bands of a Lorentzian function, respectively. We cut 9.4-9.8 /im from spectra in 
Figure [3] due to the artifacts (see text). The lower panel shows a residual spectrum at each 
plot. The bottom panel is the same as in Figure [3l 
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Fig. 4. — Continued. 
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Fig. 5. — Variation of the UIR band ratios of the 3.3 fim to the 11.3 fim band ((a) and 
(g)), the 6.2 /im to the 11.3 /im band ((b) and (h)), the 7.7 fim to the 11.3 /im band ((c) 
and (i)), the 8.6 yum to the 11.3 /im band ((d) and (j)), the 12.6yum to the 11.3 /im band 
((e) and (k)), and the 6.2yum to the 7.7 fim band ((f) and (1)) against the IRAS color of 
hbiMra/IuiMin and the 74i^74i?/ color of /l24/-^sii- The turquoise open circles and the red open 
triangles with the black errorbars indicate the UIR band ratios and the AKARI or IRAS 
colors of Groups A and B, respectively. The UIR band ratios are corrected for the extinction 
and the contribution from the hydrogen recombination lines. 
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Fig. 6. — Variation of the UIR band ratios of (a) the 7.7/im to the 11.3 fim band and (b) 
the 6.2 /im band the 11.3 fim band against the UIR band ratio of the 3.3 fim to the 11.3 fim 
band. The turquoise open circles and the red open triangles with the black errorbars indicate 
the UIR band ratios of Groups A and B, respectively. The UIR band ratios are corrected 
for the extinction and the contribution from the hydrogen recombination lines. 
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Fig. 7. — Comparison of the observed UIR band ratios with the models. The model grids 
are calculated for a mixture of neutral and ionized PAHs with a power-law size distribution 
exposed to a blackbody with different temperatures T* and various ionized fractions of PAHs 



/i. Case I (a) is calculated with (n™"^, n^^"^ 



= (20,4000) and Case II (b) with (ng'"^, n?^^^^ 



(100,4000). See Appendix for details of the model calculation. The symbols of the observed 
points are the same as in Figure O 



